AD- A 064  600  Alt  FORCE  6E0PHYSXCS  LAB  HANSCON  AFB  I 

TROPOSPHERIC  STRATOSPHERIC  TURBULENCE 
JAN  60  S  F  ZXNNERNAN'  E  A  NURFHT 
UNCLASSIFIED  AF6L-TR-60— 0010 


IA  F/6  6/1 

AND  VERTICAL  DtFFUflVXTIC—CTCIUI 


AFGL-TR-80-0020 

ENVIRONMENTAL  RESEARCH  PAPERS.  NO.  SSI 


Tropospheric  Stratospheric 
Turbulence  and  Vertical  Diffusivities 


14  January  1980 


P'l  :f.,.  r  ■  ^3'^ 

:?i;;;,;.  sJfc  frgffcg  $*3 


This  report  has  been  reviewed  by  the  ESD  Information  Office  (OI)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 


This  technical  report  has  been  reviewed  and 
is  approved  for  publication. 


FOR  THE  COMM  AN  DER 


Unclassified 


SECURITY  CL  ASS*PiC  AT(0N  0p  tHiS  PAGE  Date  Fnte  ed) 


REPORT  DOCUMENTATION  PAGE 


.RtEORT  HilMBER 


airp  READ  INSTRUCTIONS 

AU[- _ BEFORE  COS1P1.ETINT  hoRM 

1  GOVT  ACCESSION  NO.  3  RECIPIFNT’S  C  AT  A  _OC.  NUMBER 


AFGL-TR-80-0020  y  '  j/^ 

~~TTT[_Ci^J  !abiitliy  "  ""  -  _ 

T ROPOSPHERIC  STHATOSPHK  RIC 
TURBULENCE  AND  VERTICAL 
DUFF  USIVIT1ES  . 


S.  P./  Zimmerman 
E.  A.,  Murphy 

9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  ’ 

Air  Force  Geophysics  Laboratory  (LKB) 
Hanscom  AFB 
Massachusetts  01731 

1  I  CONTROLLING  OFFICE  NAME  AND  ADORES* 


5  TYPE  OP  REPORT  A  PERiOT  C  O.fBfr 

Scientific.  Final. 

6  PERFORMING  ORG  RfPC'P*  N  SvJbl R 

F  HP  No.  6!M  £ 

15  CONTRACT  op  gran’’  NjM0F«1 


10  PROGRAM  Element  pROjEC’-  Task 
AREA  A  WORK  .isriTtU  mtitos  - 


66900703 


2.  -«€Pf>WT  OAT? - 


Air  Force  Geophysics  Laboratory  (LKB)  / /  j  January  i»8t> 

Hanscom  AFB  .. _ — • 

Massachusetts  01731  IP 

T5  MONITORING  AGENCY  NAME  A  ADDRESSfrf  different  Irom  Control tin*  Office ,  15  SECuRlTV  CLASS  'of  •  .t  '**»•< 

/  K  /  Unclassified 


I  OlST  Ri 8U.XiO«-SXAT.eUfctiT  (ol  th&**p€*y  -  * 


Approved  for  public  release;  distribution  unlimited. 


i 


I  IT  DISTRIBUTION  ST  ATEMENT  (of  the  abstract  entered  fn  Block  20,  If  different  from  Report  \ 


Ms.  Supplementary  notes 


I  19.  KEY  WORDS  (Conllnu.  on  rovers*  side  If  necessary  and  identify  t>y  block  numb ei 


Turbulence 
Tropopause 
Radioactive  debris 


20  ABSTRACT  f Continue  on  reverse  elde  If  necessary  and  identity  by  block  number) 

The  determination  of  winds  and  temperature  by  rawinsonde  allows  us  to 
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Tropospheric  Stratospheric  Turbulence 
and  Vertical  Diffusivities 


I.  INTRODUCTION 


The  utilization  of  radioactive  tracers  deposited  into  the  stratosphere  to  deter¬ 
mine  effective  stratospheric -transfer  coefficients  is  well  documented  by  Reed  and 

12  3  4 

German,  Gudiksen  et  al.  Luther,  and  Danielsen  and  Louis.  In  these  studies, 

the  dispersion  of  tracers  in  the  stratosphere  and  their  rate  of  emission  from  the 

stratosphere  was  used  to  determine  "effective"  stratospheric  diffusion  coefficients. 

It  is  considered  that  in  these  two-dimensional  analyses  the  longitudinal  diffusivity 

(K  )  is  due  to  the  effect  of  mean  meridional  motions  and  large  wave  transfer.  The 

yy  6 

tilted  K  coefficients  are  due  to  the  tilt  of  the  tropopause  when  the  transfer  is 
yz 

downward,  from  the  equatorial  tropopause  to  the  polar  regions.  K  ,  the  vertical 
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diffusivity  term,  which  has  the  smallest  magnitude,  is  the  term  that  completes  the 
vertical  transfer  to  account  for  the  observations. 

The  questions  we  address  in  this  report  are:  what  is  the  function  of 
troposphere -stratospheric  turbulence  in  the  vertical  transfer  of  trace  elements, 
mass,  heat,  and  momentum,  and  how  does  it  relate  to  the  above-deduced  param¬ 
eters0  To  examine  these  questions,  we  have  utilized  a  large  bodv  of  "rawinsonde 
data"  supplied  by  the  Environmental  Data  Service  National  Climatic  Center  in 
Asheville,  North  Carolina.  This  data  consists  of  winds,  temperature,  and  pressure 
as  a  function  of  altitude  obtained  from  all  available  stations  (currently  144  stations) 
for  the  period  1948-1976.  At  the  present  time  only  the  1970-1976  data  have  been 
reformatted  for  processing  on  the  CDC  6600  computer  system.  These  twice-daily 
data  are  used  in  calculating  the  Richardson  number,  a  stability  criteria  for  deter¬ 
mining  the  presence  or  absence  of  atmospheric  turbulence,  and  is  given  by  the 
following  relationship: 

Ri  *  T  [iS  +  r]y/  [  §Z  ]  (1> 

where 

g  is  the  acceleration  of  gravity 

8T/dZ  is  the  vertical  temperature  gradient 

r  is  the  dry  adiabatic  lapse  rate  (=  9.  8  K/km) 


and 


favl  2  _  PM2  f9Vy 

I  WZ  I  -  [^Z~J 


is  the  square  of  the  vertical  shear 


of  the  horizontal  winds. 


The  rawinsonde  system  is  not  ideally  suited  for  taking  measurements  that  will 
provide  thermodynamic  data  at  the  required  sampling  intervals  for  this  application. 
The  data,  as  used  in  this  study,  were  originally  intended  to  provide  northern  hemis¬ 
phere  pressure  maps  for  a  number  of  millibar  pressure  levels.  In  order  to  deter¬ 


mine  the  Richardson  number  at  equal -height  intervals,  the  temperature  and  wind- 
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component  data  points  are  fit,  using  a  Hermite  interpolation  algorithm.  The 
temperature  function  along  with  the  derivatives  of  the  winds  and  temperature,  with 


5. 


Tsipouras,  P, ,  and  Cormier,  R.V.  (1973)  Hermite  Interpolation  Algorithm 


•pol 

for  Constructing  Reasonable  Analytic  Curves  Through  Discrete  bata  Points, 


respect  to  altitude,  are  interpolated  at  100-m  intervals  and  used  to  calculate  the 
Richardson  number.  The  probability  of  occurrence  of  turbulence  is  obtained  by 
defining  a  critical  Richardson  number  (R.  =  1/4  according  to  current  literature) 
and  assuming  that  turbulence  is  present  when  the  Richardson  number  is  equal  to 
or  less  than  this  value. 

This  approach  follows  that  of  Zimmerman  and  Murphy®  in  a  previous  article 
wherein  they  examined  the  microscale  transfer  coefficients  in  the  mesosphere 
utilizing  rocket  grenade  data.  Once  given  the  Richardson  number,  they  related 
the  ratio  of  the  vertical  turbulent  velocity  (w)  divided  by  the  mean  wind  (V)  to  this 
parameter,  following  the  empirical  relation  fitted  to  Deacon's  analvses: 


w 

V 


l  -  0.  15  (R.  >  0)  )  _ 

l  1  \  +  °- 

(  0.  15  (R.  <  0)  )  t  1 


08  . 


(2) 


Thus,  given  the  Richardson  number  and  the  horizontal  wind  velocity,  we  assume 
that  the  ratio  of  a  turbulent  parameter  to  its  mean  quantity  measured  in  the  bound¬ 
ary  layer  may  be  scaled  for  use  in  the  higher  atmosphere.  This  a  priori  assump¬ 
tion  appears  validatedby  the  resulting  amplitudes  in  the  mesosphere,  and  it  w  as  more 

recently  supported,  in  absolute  value,  by  Mansen  et  al.  Thus,  given  the  above 
2 

turbulent  intensity  <w  >  and  atmospheric  stability,  we  may  then  estimate  the 
unknown  turbulent  parameters  described  below. 


2.  STRATOSPHERIC  AND  MESOSPHERIC  TURBULENCE 

The  assumptions  used  in  the  following  analysis  is  that  once  we  are  given  the 
2 

turbulent  intensities  <w  >  from  Eq.  (2)  we  can  assume  that  the  spectrum  of  the 
turbulent  motions  will  be  inertial  and,  as  discussed  below,  the  vertical  length  scale 
will  be  limited  by  buoyance.  Thus,  given  inertiality  and  one  dimensionality,  we  have 

OC 

<Vb2>  =  f  E  (k)  dk  (3) 

and 


6.  Zimmerman,  S.P.,  and  Murphy,  E.A.  ( 1977)  Stratosphere  and  mesospheric 

turbulence  in  Dynamical  and  Chemical  Coupling,  D.  Reidel,  Dordrecht, 
Holland,  pp  35-47. 

7.  Deacon,  E.  L.  (1959)  The  problem  of  atmospheric  diffusion,  Int.  J.  Air  Pollut., 

pp  92-  108. 

8.  Manson,  A.  H.  (1980)  to  be  published. 


, .  , .  ,  2  /  3  ,  -  5  /  3 

h  (ki  =  at  k 

2  2 

where  <v  >  =  <  w  >  is  the  buovancy-limited,  vertical  turbulent  intensit\  -et 
“  2 

identically  equal  to  <w  >,  as  given  in  Eq.  (2);  K  < k )  is  the  turbulent  energv 
spectrum;  a  is  the  Kolmagoroff  constant,  and  for  the  vertical  dimension  it  has  the 
value  '  0.  8;  k  is  the  wavenumber;  and  e  is  the  rate  of  turbulent  dissipation.  Thus, 
from  Eq.  (3)  we  derive 


„  2/3  -2/3 

kB 


(4) 


In  a  review  of  the  Lumley-Shur  theory  of  turbulence  in  a  buoyant  stably  stratified 
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medium,  it  is  stated  that  the  turbulent  spectrum  will  be  inertial  in  the  high  wave- 
number  region,  and  that  the  "buoyant"  subrange  will  strictly  be  a  function  of  at¬ 
mospheric  stability  in  the  form  of  the  Brunt  Vaisalla  frequency  (N)  and  wavenumber 
(k).  As  a  consequence  of  this,  t ,  which  is  conserved  in  the  inertial  region,  is  a 
function  of  wavenumber  in  the  buoyant  region.  Thus,  the  spectral  energy  is  given 
by 


E  (k)  =  Cj  N2  k~3  +  at  2/3  k'5/3 


=  at 


2/3 


-5/3 


(5) 


where  k_  is  the  transition  wavenumber 


between  the  inertial  and  buoyant  subranges.  In  this  report,  we  presume  that  the 
vertical  turbulent  spectrum  is  "inertial"  and  limited  in  spectral  extent  to  the  tran¬ 
sition  wavenumber  (kg),  and,  thus,  the  vertical  turbulent  kinetic  energy  as  given 
by  Eq.  (4)  becomes 


9.  Phillips,  O.  M.  (1967)  On  the  Bolgiano  and  Lumley-Shur  theories  of  the 

buoyancy  subrange,  in  Atmospheric  Turbulence  and  Radio  Wave  Propagation, 
Nauka,  Moscow,  USSR,  pp  121-128. 
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For  a  =  ~  0.  8  and  Cj  -»  2  to  2.  5  (for  simplicity  say  2.  5),  we  have 


2 

<  w  > 


1  t 

7  N 


(7) 


or 


e  =  2  <  w  >  N 


Heisenberg111  had  proposed  an  effective  diffusivity  (K)  in  explaining  the  establish¬ 
ment  of  an  inertial  subrange  (Hinze1  )  as 


oo 


1/2 

dk  . 


(8) 


This  reduces  to 


~  1  £  1  <w2> 

B  ‘  4  ^2T  =  7  “TT-  ■ 


(9) 


It  is  interesting  to  note  here  that  there  is  a  factor  of  approximately  1/2  intro¬ 
duced  into  the  amplitude  of  the  turbulent  parameter  £q  via  the  definition  of  the 
turbulent  spectrum  Eq.  (5).  If  instead  we  had  used  the  rationale  that  the  one- 
dimensional  spectrum  is  only  considered  to  the  buoyant  limit,  then,  Eq.  (5)  would 
read 


EM,at2l3tB213 


and  Eqs.  (7)  and  (9)  would  follow  as 


✓  2  -^c  ~  1  <  w  > 

<w  >.  wandKB=  Ir¬ 


respectively. 


X  DATA  AND  RESULTS 

In  this  preliminary  report,  we  examine  the  seasonal  dependence  of  the  vertical 
turbulent  diffusivity  (K)  as  a  function  of  latitude  and  altitude,  covering  the  latitude 
band  12°Nto7I®N.  In  a  future  report,  we  shall  present  in  the  same  manner: 


10.  Heisenberg,  W.  (1948)  Z.  Physik,  pp  124-628. 

11.  Hinze,  J.  O.  (1959)  Turbulence,  McGraw-Hill,  New  York. 
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(1)  occurrences  for  R  -  0.  25.  (2)  the  resultant  turbulent  intensities  <w“'»  .  anil 
(3)  the  rates  of  dissipation  of  turbulent  kinetic  energy  u  )  over  a  larger  latitudinal 
region. 

Figure  la  is  a  representative  set  of  rawinsonde  east-west  wind  component  data. 
The  symbols  are  the  discrete  data  values  and  the  solid  curve  is  t  he  output  of  the  Hermite 
interpolation  algorithm  developed  bv  Tsipouras  and  Cormier,  3  and  here  simply  re¬ 
ferred  to  as  a  "cubic  spline  fit".  In  Figure  lb  the  average  shear  taken  from  a  linear  fit  to 
the  data  values  is  compared  to  the  derivative  taken  from  the  Hermite  interpolation 
algorithm,  that  is,  using  the  derivative  of  the  polvnominal  with  the  coefficients  of 
the  fit.  Using  this  technique,  a  close  approximation  to  the  average  wind-component 
gradients,  or  shears,  and  the  temperature  gradients  are  obtained  at  regular  inter¬ 
vals,  by  computer,  without  introducing  problems  associated  with  extremely  large 
or  small  numbers.  The  gradient  Richardson  number  determined  in  this  wav,  using 
the  average  gradients,  provides  only  a  relative  measure  of  atmospheric  stability 
or  instability. 


WIND  SPEED  (m/ttcl 

Figure  la.  Cubic  Spline  Fit  of  the  Rawinsonde  East-West  Wind  Data.  The  solid 
line  is  the  splined  fit  and  the  crosses  (X)  are  the  data  points 
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Figure  lb.  The  East-West  Wind  Shear.  The  linear  calculations  and  the  cubic 
spline  fit  are  discussed  in  the  text 


The  database  covers  the  years  1972  to  197  6,  and  for  this  preliminary  report 
we  shall  basically  consider  the  analysis  of  the  year  197  6.  The  locations  chosen  are 
Barrow,  Alaska  (Lat  7  1  °  1 8 '  N),  Yakutat,  Alaska  (59°31'  N),  Dayton,  Ohio  (39°52'  N), 
Waycross,  Georgia  (3 1  °  15'  N),  Grand  Cayman,  B.W.l.  (19°18'  N),  and  San  Andres, 
Columbia  ( 1 2 0 3 5 '  N).  Figures  2a,  2b,  and  2c  demonstrate  the  altitude  distribution 
of  turbulent  diffusivities  [Eq.  (9)]  for  three  of  these  sites  for  the  spring  and  winter 
seasons  from  1  km  to  ~  30  km  altitude.  The  data  are  calculated  at  100-m  intervals 
from  the  interpolated  wind  and  temperature  profiles  and  then  averaged  over  1-km 
intervals.  What  is  immediately  obvious  is  the  striking  reduction  of  turbulence  at 
the  tropopause  (noted  on  the  figures  by  the  solid  dots)  and  the  significant  enhance¬ 
ment  of  winter  diffusivities  over  those  of  the  spring  season,  particularly  so  in  the 
stratosphere. 

As  observed  for  those  latitudes  when  the  tropopause  is  unambiguously  deter¬ 
mined  (Figure  2a,  San  Andres  12°35'  N),  the  rapid  reduction  of  the  turbulent  diffu¬ 
sion  coefficient  at  the  tropopause  altitude  clearly  delineates  the  disturbed  region  of 
the  lower  atmosphere  from  the  more  placid  stratosphere.  A  similar  phenomenon 
is  observed  at  Barrow  (70°  18'  N),  but  the  tropopause  is  difficult  to  determine  using 
the  conventional  logic.  At  Great  Falls,  there  is  again  a  marked  reduction  of  tur¬ 
bulence  at  the  tropopause,  but  significant  turbulence  in  the  stratosphere  for  the 


winter  season.  The  spring  season,  on  the  other  hand,  at  all  sites  shows  very  little 
turbulence  in  the  stratosphere,  even  at  the  midlatutude  sites. 

The  rapid  reduction  of  turbulent  diffusion  is  an  unambiguous  observation  that 
suggests  a  more  appropriate  indicator  of  the  tropopause  would  be  the  layer  where 
there  appears,  in  essence,  a  step  function  of  a  parameter  involving  both  dynamical 
and  thermal  structure  rather  than  a  thermal  demarcation  alone.  The  turbulent 
diffusion  coefficient,  as  calculated  here,  is  such  a  parameter,  clearly  delineating 
a  relatively  placid  region  from  the  more  disturbed  region.  The  strong  corres¬ 
pondence  between  the  tropopause,  when  determinable,  and  the  delineation  altitude 
determined  by  the  rapid  cessation  of  turbulence  is  clearly  displayed  in  Figures  3a 

and  3b.  They  are  also  contrasted  to  the  isodensity  contours  of  radioactive  debris 

4 

as  reported  by  Danielson  and  Louis,  and  it  is  clear  that  the  tropopause,  deter¬ 
mined  by  the  cessation  of  turbulence,  shows  a  remarkable  correspondence  to  that 
level  of  the  atmosphere  where  the  trace  particles  are  rapidly  depleted. 


Figure  2a.  One  Kilometer  Average  of  the  Vertical 
Turbulent  Diffusivlties  (Kzz>  as  Given  by  Equation  (9) 
for  the  1976  Winter  and  Spring  Seasons  at  San  Andres, 
Columbia 
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Figure  2b.  One 
Kilometer  Average  of 
the  Vertical  Turbulent 
Diffusivities  (KZ7)  as 
Given  by  Equation  (9) 
for  the  1976  Winter 
and  Spring  Seasons  at 
Great  Falls,  Montana 


Figure  2c.  One 
Kilometer  Average  of 
the  Vertical  Turbulent 
Diffusivities  (Kzz)  as 
Given  by  Equation  (9) 
for  the  1976  Winter 
and  Spring  Seasons  at 
Barrow,  Alaska 


3 


90 


SO 


0 


60 


90 


60 

NORTH 


OBSERVE! 


J__l _ 1 _ I _ L- 

so 

south 


Figure  3a.  Comparison  of  the  Isodensity  Contours 
of  Radioactive  Debris  (heavy  solid  lines)  With  the 
Isocontours  of  Vertical  Turbulent  Diffusivities  in 
the  Northern  Hemisphere  (light  solid  and  dashed 
lines).  The  numbers  on  the  contours  of  turbulent 
diffusivity  are  values  expressed  in  K  X  10“^ 
cm^ /sec 


Figure  3b.  Comparison  of  the  Isodensity  Contours 
With  Three  Parameters  Marking  the  Tropopause 


The  large  variability  of  stratospheric  turbulence  is  further  demonstrated  in 
Figure  4  where  we  observe  a  fairly  intense  increase  of  turbulence  in  the  stratos¬ 
phere  for  the  winter  season  of  1976,  as  compared  to  that  for  1975.  The  tropos¬ 
pheric  turbulence  in  1975,  on  the  other  hand,  is  only  slightly  enhanced  over  that  of 
1976.  This  large  variability  of  the  stratospheric  turbulent  vertical  diffusion  above 
the  tropopause  will  have  a  significant  effect  upon  seasonal  and  yearly  variability  of 
transfer  across  the  tropopause. 
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GREAT  FALLS,  MONTANA 
47*29'  1975 
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TURBULENT  DIFFUStVITY  (K)(m?/sec) 

Figure  4.  The  Vertical  Turbulent  Diffusivities  of 
the  1975  Winter  and  Spring  Seasons.  These  are  for 
comparison  with  the  1976  results.  Figure  2c,  to 
demonstrate  the  large  yearly  differences  observed 
in  the  stratosphere 


Effective  vertical  transfer  has  been  derived  by  the  diffusion  analysis  of  the 

1-3 

stratospheric  decay  and  tropospheric  deposition  of  radioactive  bomb  debris. 

1 2 

Here  the  meridional  transient  eddies  generated  by  Murakami  were  utilized 
(References  1-3)  to  determine  the  meridional  large-scale  diffusivities,  and  then  the 
continuity  equation  was  used  to  estimate  a  vertical  diffusivity  that  best  approxi- 

4 

mated  the  observations  cited  by  Danielson  and  Louis.  The  comparison  of  their 
results  with  the  measurements  is  given  in  Figure  5  (their  Figure  9.  11),  and  it  is 
quite  obvious  that  the  derived  slopes  of  the  isodensity  tracer  contours  are  not  at 
all  commensurate  with  the  measurements  in  the  northern  latitude.  This  conflict 


with  the  observations  could  be  due  to  a  number  of  facets,  one  being  incorrect 
values  of  KyZ>  the  diagonal  element  of  the  meridional  transfer  function  based  upon 
the  large-scale,  heat -transfer  calculation.  Another  is  in  their  derived  vertical 


diffusivity.  This  disagreement  is  evident  in  the  comparison  of  our  Figure  3a  and 


12.  Murakami,  T.  (1962)  Stratospheric  Wind  Temperature  and  Isobaric  Height 
Conditions  During  the  1.  (j.  Y.  Period,  Parti,  M.  I.  T.  Report  No.  5 
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their  Table  9.  1  (reproduced  here  as  Figure  6)  that  shows  the  iso  turbulent  diffusion 

4  2  4  2 

coefficients  lines  for  K  =  1  X  10  cm  1  sec  and  =  1.  5  >  10  cm  /sec.  This  later 

isocontour  is  also  compared  to  the  observations.  There  is  no  parallelism  to  the 
3 

10  PC./m  curve  as  this  analysis  displays,  but  a  deeper  penetration  into  the  mid¬ 
latitude  stratosphere  with  larger  values  than  we  determine.  Our  tropospheric 
values,  on  the  other  hand,  are  somewhat  larger  than  theirs,  approaching  an  order 
of  magnitude  difference  when  the  similar  averaging  distance  (3  km)  is  approxi¬ 
mated  (Figure  3a). 


Figure  5.  Observed  Radioactivity  Isodensity  Contours 
(dashed  lines)  as  Compared  to  the  Model  of  Danielsen 
and  Louis 


Figure  6.  Spring  Values  of  the  Isocontours  of  the 
Vertical  Turbulent  Diffusivities  (light  solid  and  dashed 
lines)  as  Contrasted  to  the  Observed  Radioactivity 
Isodensity  Contours 
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4.  CONFUSION 

In  conclusion,  we  have  demonstrated  that  vertical  turbulent  transfer,  as 
determined  from  rawinsonde  Richardson  number  calculations  and  utilization  of 
boundary  layer  analysis,  results  in  average  turbulent  diffusion  coefficients  in  the 
tropopause  and  lower  stratosphere.  Comparison  of  iso-diffusion  contours  with 
radioactive  debris  isodensity  contours  demonstrates  remarkable  parallelism. 

Since  turbulent  diffusion  operates  upon  species  gradients  (that  is,  v^.jj.  «  K  (-^  ~ )), 
we  may  infer  that,  indeed,  the  experimentally  determined  diffusivities  dictate  the 
maintenance  of  a  flux  perpendicular  to  the  species  isocontours  and,  thus,  continued 
maintenance  of  these  seasonal  contours.  Future  efforts  will  be  to  extrapolate  the 
results  to  southern  latitudes  and  demonstrate  other  seasonal  effects,  and  continue 
comparisons  with  other  data. 
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